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ABSTRACT 

We use a series of idealized, numerical SPH simulations to study the formation and evolution 
of galactic, gas-rich disks forming from gas infall within dark matter halos. The temperature and 
density structure of the gas is varied in order to differentiate between (i) simultaneous gas infall at 
a large range of radii and (ii) the inside-out build-up of a disk. In all cases, the disks go through 
phases of ring formation, gravitational instability and break-up into massive clumps. Ring formation 
can be enhanced by a focal point effect. The position of the ring is determined by the angular 
momentum distribution of the material it forms from. We study the ring and clump morphologies, 
the characteristic properties of the resulting velocity dispersion field and the effect of star formation. 
In the early phases, gas accretion leads to a high vertical velocity dispersion. We find that the 
disk fragmentation by gravitational instability and the subsequent clump-clump interactions drive 
high velocity dispersions mainly in the plane of the disk while at the same time the vertical velocity 
dispersion dissipates. The result is a strong variation of the line-of-sight velocity dispersion with 
inclination angle. For a face-on view, clumps appear as minima in the (vertical) dispersion, whereas 
for a more edge-on view, they tend to correspond to maxima. There exists observational evidence of 
a systematic variation of the velocity dispersion with inclination angle in high-redshift disks, which 
could be partly explained by our simulation results. Additional energetic sources to drive velocity 
dispersion that are not included in our models are also expected to contribute to the observational 
results. 

Subject headings: galaxies: high-redshift — galaxies: evolution — galaxies: irregular — galaxies: 
formation — galaxies: starburst 



1. INTRODUCTION 

Recent observations at redshift z ^ 2, close to the peak 
of the cosmic star formation history, have revealed a sig- 
nificant population of massive (baryonic mass Mbar '^ 
lO^^M©), extended [R ^ Q kpc) disk galaxies with star 
forma tion rates of the order of SFR ^ 100 M© yr~^ 
(e.g. iForster Schreiber et al.l l2006l [20091 ). The corre- 
sponding observed gas velocity dispersions are very high, 
cr ~ 40 — 90 km s~^ and the rotation-t o-dispersion ratios 
y/a ~ 2 - 6 thus r elativ ely low (e.g. iCresci et al.1l2009l: 
iGenzel et al.l 12006'. '2008") . In addition, the galax ies are 
gas- r i ch with typical gas fractions of ~ 50% (Dadd i~et al.l 
l2010t iTacconi et al.ll20l0[ ). Star formation within them 
takes place in a few giant clumps with typical masses 
of Mniiimp ^ 10^M(T) and radii of i?ciump '^ 1 kpc 
(jElmegreen et al.l |2009|) . In some of the sources rings 
have been detected, suggesting an evolutionary sequence 
from rings to clumps to bulges on a timescale of t < 1 Gyr 
(JGenzel et al. 20081 . 

The properties of this observed z ~ 2 galaxy pop- 
ulation are very different from present-day disk galax- 
ies. Typical disk galaxies, such as the Milky Way, dis- 
play smooth surface density profiles with no evidence 
for supermassive clumps. The difference to the z ~ 2 
population can be illustrated by the typical parame- 
ters for the local disk popul ation of SFR '^ 1 M© yr~^ 
(jRobitaille fc Whitnevll2010l and references therein) and 



V/a - 20 (iDib et al.l 120061) . Star formation takes place 
in giant molecular clouds with typical masses of Mgmc ~ 
lO^Af0 and has been constant within a fa ctor of a few 
for the last t ~ 10 Gyr {.Aumer fc Binnev 2009). As- 
sociations of molecular clouds can have masses as high 
as a few lO^M© (e.g. iRand fc Kui^ml Il990l ). which 
is the largest scale of fragmentation in z ~ galaxies. 
However, the thick disk components of low-redshift disks 
exhibit stellar ages consistent with a formation at z ^ 2 
and high velocity dispersions (e.g. iSchonrich fc BinnevI 
|2009)- It has been proposed that the thick disks have 
forme d in massive, clumpy high-redshift disks with high 
SFRs (IForster Schreiber etal.l [200l IGenzel et a"Ll l2006l: 
iBournaud et al.ll2007l I2OO90 . 

Observations indicate that only a smaller fraction (~ 
1/3) of the z ~ 2 galaxies are obvious major merger s 
(jForster Schreiber et all2006l[2009tlShapiro et al.l2 008'). 
Thus, the high star formation rates and high velocity 
dispersions seem not to be merger-induced in the ma- 
jority of the sources. Instead, gravitational instability 
of differentially rotating disks with high surface densities 
(jToomrd 11964) has been identified as a key process for 
the fragmentation of z ~ 2 disk galaxies, which are then 
likely to evolve into bulge -dominat ed svstems (iNoguchil 
[19991: llmmeh et al.ll2004al ). Lately. lRo5leo et al.l (|20iol ) 
extended the classical Toomre instability criterion to tur- 
bulent disks which show a rich variety of global instabili- 
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ties, providing a theoretical framework for gas-rich high-z 
disk evolution. 

Clumpy and turbulent disks resembling the z ~ 2 
galaxies have recently also been produc ed in cosmo- 
logical simulations of galax y formation (jAgertz et alJ 
I2009bt iCeverino etUI l2010[ ). As argued by the au- 
thors, key aspects in the success of these simulations 
are a high spatial resolution and a cooling description 
that follows the thermal evolution of the gas below 
T < 10^ K allowing a multi-phase ISM to develop. 
In these simulations, the high gas fractions and star- 
formation rates are power ed by cold streams of gas mov- 
ing along the cosmic web (iKeres et al.ll2005l : iDekel et al.l 
I2009bt I Johansson et al.l l2009bD . However, the physical 
source of the large observed irregular motions is still un- 
clear. Several scenarios have be en proposed in the liter- 
ature including acc r etion (e.g. Forster Schrcibe r et al.l 
20061 iGenzel et al.l [2001 [Elmegreen fc Burkerti |201( 
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Klessen fc Hennebehel |2010D. stellar radiative feedback 
(jMurrav et al.l 12010!: 'Krumholz fc Dekell |2010D and disk 



self-gravity (e.g. Immcli ct al. 2004a; Bournaud et al.l 
[2007[ iBurkert et alJl201Q) . 

In this paper we present SPH simulations of idealized 
disk formation by gas infall via a cooling flow within a 
galactic halo. We analyze the evolution of the morphol- 
ogy and the turbulent structures in the forming disks. 
Both models with and without star formation and SN 
feedback and with different initial temperature and den- 
sity profiles are presented. We demonstrate that ring for- 
mation and subsequent gravitational fragmentation oc- 
curs in all of the scenarios. In addition, we analyze the 
structure of velocity dispersion created by accretion and 
disk self-gravity and finally discuss the efficiency of stel- 
lar disk heating in our simulations. 

We describe the setup of our simulations in Section 
2, present our results in Section 3 and discuss them in 
Section 4. 

2. SIMULATIONS 

The simulations were perform ed using the TreeSPH- 
code GADGET-2 (|Springell l2005') on the local Ahix 3700 
Bx2 machine. The code includes the standard radiative 
cooling rates for an optically thin, primordial composi- 
tion of h ydrogen and heliu m in ionization equilibrium 
following iKatz et al.l (|1996l ). Their figure 1 displays the 
cooling rate as a function of temperature. The cooling 
curve features a cutoff below T ~ lO^K and two peaks 
at T - lO-^-^K and T - 10^"K resuhing from colhsional 
excitation of H and He, respectively. We included a spa- 
tially uniform tim e-independent UV backg round, appro- 
priate for z = 2 (jHaardt fc Madaul Il996f ). The abun- 
dances of the different ionic H and He species were com- 
puted by solving the network of equilibrium equations 
self-consistently for a specified value of the UV back- 
ground radiation field. 

The chemical composition of the infalling gas could 
possibly be pre-enriched. However, the actual metal con- 
tent is difficult to estimate, as there are no direct obser- 
vations of infalling gas at 2 ~ 2 (e.g. Steidcl ct al. 201(J 
and observed galaxies and outfiows are likely to be signif- 
icantly more metal-rich than the infalling gas. Consider- 
ing the metallicity of t yoical damped Lym an a systems 
at z ~ 2, Z ~ 0.1^0 (|Pettini et al.lll994[ ). the cooling 
rates would be enhanced for temperatures T > W^-^K 



by a factor of a few ([Sutherland fc Dopitalll993l ). In our 
simulations, this would lead to an increase in gas mass 
flow to the center of the halo. However, due to the large 
uncertainties in the actual gas metallicity at 2 ^ 2, we 
rely on the assumption of a primordial chemical compo- 
sition in the present study. 

Star formation and the associated supernova feedback, 
when included, are m odeled using the prescriptions of 
iSpringel fc Hernguistl (jJOOS). This sub- resolution model 
uses a statistical formulation to deal with physics act- 
ing on scales, which are not resolved by simulations of 
galaxy formation. Above a threshold density pth , each 
SPH particle is assumed to represent a fluid comprised 
of cold, condensed clouds in pressure equilibrium with 
an ambient hot gas. Such a multiphase treat ment of the 
interst ellar me dium is motivated by McKcc fc Ostriken 
(J1977!) (see also [Johansson fc Efstathioul l2006'). A steUar 
population is modeled to form from the cold medium on 
a timescale ts , which is proportional to the dynamical 
timescale of the gas, so that tg = is,o(p/pth)~'^'^ • The 
energy released by one SNII is typically 10^^ ergs per su- 
pernova, which may be expressed as a supernova temper- 
ature TsN — 2/iUsN/(3fcB) ~ lO^K. The hot phase is re- 
quired to have a temperature in excess of lO^K, whereas 
the cold phase is assumed to be at Tc ~ lO^K. The 
threshold density, pth , is determined self-consistently in 
the model by requiring that the equation-of-state (EOS) 
is continuous at the onset of star formation. Stars form 
from the cold clouds in regions were the number den- 
sity n > nth = 0.128cm~'^. The parameters govern- 
ing the mode l are set to repro duce the observed Kenni- 
cutt relation (|Kennicuttlll998[ ) , yielding a star formation 
timescale tgfi = 2.1Gyr. The unphysically low star for- 
mation density threshold density nth is a consequence of 
the relatively low resolution of typical galaxy formation 
simulations, such as ours (~ 10^ Mc^ per partic le). As 
has recently been shown ([Governato et al.ll2010( ). higher 
density thresholds are appropriate for higher-resolution 
simulations of dwarf galaxies (~ 10'^ M0 per particle). 

Radiation pressure, which might play an important 
role in high-2 star forming galaxies (see iMurrav et al.l 
I2010D . is not included in our simulations. 

The setup fo r our simu lations was motivated by 
iKaufmann et al.l ([20061 120071 ). We study the evolution of 
an initially hot baryonic component within a dark matter 
halo with an NEW density profile 
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([Navarro et al.|[r996D . The str uctural parameters of the 
halo were chosen according to iZhao et al.l ([2OO90 . The 
virial radius i?vir and virial mass Mvir are defined ac- 
cording to the spherical virialization criterion. The cor- 
responding virial over-density Avir(2 ~ 2) « 187 de- 
pends on redshift 2 and on the assumed cosmology. We 
adopt M yir = lO^^Mfyi as exp ected for 2 '-^ 2 star-forming 
galaxies ([Dekel et al.ll2009b[ ). Eor this halo mass, 2 = 2 
and a ACDM cosmolo gy with par ameters according to 
iHinshawetaD ([20091 ) . IZhao et al.l ([2009( ) find a concen- 
tration parameter of c = R^-^/Rs = 4.2 and a virial 
radius of i?vir = 108 kpc, which we adopt. We apply a 
fiducial virial spin parameter of Adm = /|°" = 0.08 
([Bullock et al.l I2001D . where jdm is the average spe- 
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cific angular momentum of the dark matter within i?vir 
and V = y/GMviv/Rvii- is the halo circular velocity at 
i?vir- This value for Adm is on the high end of the dis- 
tribution of spin pa rameters of CDM halos found by 
I Bullock et alJ (|2001| ). However, similarly high spin pa- 
ramet ers have been inferre d for extended disk galaxies at 
z ~ 2 (jBurkert et al.l[20Tol ). thus motivating our choice. 

The halo is set up using a population of TVdm = 
10^ collisionle s s pa rticles follo wing the description o f 
iSpringel et alJ (2005) (see also [Johansson et al.ll2009al ). 
We truncate the halo at 1.5 i?vir. The particle mass is 
todm = 1.34 X 1O^M0 and the gravitational softening 
length is set to edm = 500 pc. 

A gas component was included with a baryonic-to-dark 
mass ratio of 15%. The gas component is also truncated 
at 1.5 i?vir • We use iVbar — 10^ particles, a particle 
mass TTibar = 2.01 X 1O^M0 and a gravitational softening 
length of Cbar = 250 pc . The dominant mode of gas 
accre tion onto forming g alaxies is a topic still under de- 
bate (iDekel et al.ll2009al) . Smooth inflows of cooling gas, 
that has previously been shock-heated, cold streams and 
mergers are among the discussed modes and are all likely 
to play a role. We intend to keep our models simple in 
this context and therefore assume spherically symmetric 
profiles of density pbar and temperature T. The varia- 
tion of these profiles allows us to study different scenarios 
of infall. In our models, the gas component is initially 
hot (T > lO^K), with an internal energy structure deter- 
mined by the assumpt ion of hydrostatic equilibrium (cf. 
iKaufmann et al.ll2006( ). The pressure profile is thus 
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Unlike Kaufmann et al., we do not impose a lower tem- 
perature floor, as we do not want to prevent the forming 
disk from fragmenting. The build-up of the disk is de- 
termined by the initial profiles of gas density, internal 
energy and angular momentum. We vary the initial con- 
ditions as follows: 

• Model A assumes that the gas component is ini- 
tially well mixed with the dark matter component. 
Thus, the gas density also follows an NFW density 
profile similarly to the DM component, but with a 
correspondingly lower density normalization ps (cf. 
IKaufmann et al.| [2006'). The gas is in hydrostatic 
equilibrium with an initial temperature of the order 
of T = 10^ K within i? < 20 kpc. The initial spe- 
cific angular momentu m profil e is a power la w with 
7 PC R (cf. [Kaufma nn et al.l [2006: Bulloc k~et al.l 

|2001[ ) , i.e. the rotational velocity is independent 
of the radius. 

• With Model B we intend to explore significantly 
different central initial conditions and to conse- 
quently alter the initial formation stages. We as- 
sume a constant gas density within i? < 20 kpc, 
adjusted so that the total mass within this radius is 
equal to that of Model A, Mfl,<2o = 2.5 x W^Mq. 
The temperature in this model at _R = 20 kpc is 
also of the order of T = 10^ K. However, due to 
the assumption of hydrostatic equilibrium the tem- 
perature increases with decreasing radius resulting 
in a factor of ~ 5 higher temperature in the center. 
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Fig. 1. — Cooling timescales tcool a-s a function of radius r for 
Models A (blue), B (red) and C (green). The solid lines are for the 
initial conditions, whereas the dotted lines are for simulation time 
t = 200 Myr. 



thus creating a 'hot bubble'. In the constant den- 
sity region the angular velocity of the gas is set to 
a constant value. Outside i? = 30 kpc the density 
and angular momentum profiles are as in Model A 
with a smooth transition region between 20 and 30 
kpc. 

• Model C has identical gas density and angular mo- 
mentum profiles to Model B, but the condition 
of hydrostatic equilibrium is removed in order to 
model an accretion history that is intermediate to 
Models A and B. Within i? < 15 kpc the tempera- 
ture is set to a constant value of T ~ 3 x 10^ K, be- 
ing lower than the central temperature of Model B, 
but higher than the central temperature of Model 
A. Between 15 and 30 kpc the temperature de- 
creases linearly with radius and outside i? > 30 kpc 
the temperature profile is identical to Model B. 

• With Model SF we intend to study the ef- 
fects of star formation and supernova feedback 
on the ev olution of ou r system s. As has been 
shown by llmmeli et al.l (|2004a|) . the specific im- 
plementation of these processes has a large influ- 
ence on the formation and evolution of unstable 
disks. We therefore rely on the generic model of 
lSpringel"fc Hernquist (2003J) to be able to draw 
general conclusions. We choose identical initial 
conditions to Model B in order to make a direct 
comparison possible. 
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"g^" ~ '7Wr~' ^^^'^'^ Jgas is the 
average specific angular momentum of the gas within _Rvir 
is for all models set to a fiducial value of 0.08, so that 
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Fig. 2. — A cut perpendicular to the disk region (y=0) of Model 
A at i = 35 Myr. The color represents density p and over-plotted 
in yellow is the velocity field illustrating the focusing effect leading 
to accretion into the ring region. The velocity increases from a fewr 
10 kms"^ to ~ 160 kms~^ near the disk. 



3. RESULTS 

Due to cooling the baryonic component evolves out 
of hydrostatic equilibrium rapidly and begins to form a 
disk in the center of the halo. The cooling timescale, 
icooi oc p~^ and the temperature, density and angular 
momentum profiles determine the evolution and struc- 
ture of the forming baryonic component. In Fig. [T] we 
plot the cooling timescale as a function of radius R for the 
three different initial conditions of Models A/B/C and 
for the evolved systems at simulation time t — 200 Myr 
in order to illustrate the significant differences between 
the three models. 

As already noted in iKaufmann et all (|2006[ ) thermal 
instabilities (jBurkert fc LinI l2000t ) caused by numerical 
noise in the initial conditions lead to the formation of 
cold, dense clouds in the halo which fall onto the disk 
through a hotter medium. This however does not af- 
fect the fragmentation of our disks into clumps. The 
clouds typically consist of a number of SPH particles 
A'cioud, which is smaller than the number of particles 
in the smoothing kernel iVgpH = 40 and are thus not 
resolved. 

3.1. How the Models differ 

Our different initial conditions have been constructed 
in order to study different early phases of disk growth. 

• In Model A the cooling time increases monotoni- 
cally as a function of increasing radius. The already 
relatively dense central gas (i? ^ 4 kpc) cools very 



rapidly on a timescale of ic 



< 



20 Myr (cf. Fig. 



[T]) and forms a rotationally stabilized inner disk. 



Material within i? < 1 kpc of the rotation axis is 
in centrifugal equilibrium and falls vertically onto 
the disk, whereas the remaining gas falls in radially 
(cf Fig. [2]). Because of the peaked density profile 
and correspondingly short collapse timescales, the 
accretion of material that was originally close to 
the rotation axis and therefore in centrifugal equi- 
librium soon becomes unimportant. Now gas with 
higher angular momentum from the outer regions 
begins to settle into the outer regions of the disk, 
leading to an inside-out growth of the disk. In gen- 
eral, in this model, accretion predominantly takes 
place in the outer disk region. 

• According to equation [2] the temperature decreases 
with increasing radius in the constant density re- 
gion of Model B, thus the outer regions of the con- 
stant density sphere have the shortest cooling times 
^cooi ~ lOOMyr (cf. Fig. [T]) and the first dense 
and cold structures are formed here. These clouds 
are gravitationally accelerated inwards, moving to- 
wards the center through the hotter inner gas lay- 
ers, which subsequently also cool and undergo the 
same process. The result is that almost all of the 
material from the initial constant density region 
settles in a disk at a relatively large range of radii 
(_R < 5 kpc) within a relatively small time interval 
(t - 200 - 260 Myr). The formation phase of the 
central disk is thus distinctly different from Model 
A. 

• In Model SF star formation sets in after the first 
denser structures have formed due to cooling and 
before they enter the disk region. However star 
formation rates in this phase are low compared to 
the gas accretion, leading to an increase in the 
disk's gas mass. The inclusion of feedback in- 
creases the effective cooling timescales as the type 
II supernovae add energy to the star-forming gas. 
The central build-up of the disk is thus slightly 
delayed (i > 220 Myr) but otherwise similar to 
Model B. The first row of Fig. [3] displays the face- 
and edge-on surface density maps of Model SF at 
t ~ 230 Myr. 

• Because of the constant density and temperature 
of the gas within the central R < 15 kpc in 
Model C the cooling time is constant everywhere, 
^cooi ^ 200Myr (cf. Fig. [T]). After a first quasi- 
static cooling phase, infall starts basically every- 
where at the same time. The free-fall timescale in- 
creases with radius, resulting in an inside-out disk 
formation starting at i ~ 140 Myr. Compared to 
Model A the different density and angular momen- 
tum profiles result in the accretion into the cen- 
tral regions being significant over a much longer 
timescale (until t ^ 320 Myr). 

The dashed lines in Fig. [T] show the cooling time pro- 
files for simulation time t — 200 Myr for the Models 
A/B/C. For Models A and B, extended central struc- 
tures have formed, indicated by the short cooling times. 
The material is effectively at the lower T limit of the cool- 
ing function at T ~ W^K. Outside 10 kpc the cooling 
rates are now larger than initially, as the initial material 
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Fig. 3. — The variety of morphologies displayed by Model SF. Prom left to right: Surface density S maps of: gas face-on, gas edge-on, 
stars face-on and stars edge-on. From top to bottom: Model SF at 230, 300, 370, 450, 1000, 1500 Myr. 
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Fig. 4. — From left to right: surface density S vs disk radius R, the azimuthally averaged gravitational potential $ vs disk radius R and 
a surface density map of the disk regions. The vertical blue lines in the left and central panels mark the position of the ring. From top to 
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at these positions has moved inwards and hotter mate- 
rial from the outskirts is falling in. For Model C, at 200 
Myr, the central object is still relatively small, which is 
why its cooling time profile differs from those of Models 
A and B at this time. 

3.2. Ring formation 

All of our models exhibit ring structures and a phase 
dominated by massive clumps. The ring is very clearly 
seen in Models A and B, but is significantly less distinct 
and more obscure in Model C. The ring formation radius 
varies from model to model as displayed in the left panels 
of Fig. 12 

• In Model A a shallow peak in the radial surface den- 
sity distribution is produced soon after the start of 
the simulation at the equilibrium position between 
the halo gravity and the initial angular momen- 
tum of the gas contributing to the first stage of 
disk formation (t< 100 Myr). New disk material 
during this phase preferentially settles at the outer 
region of the already existing disk. Thus, the edge 
of the disk is always relatively distinct, resulting 
in a steep rise of the gravitational potential in this 
region. This in turn focuses the accretion of ma- 
terial to this region (cf. Fig. [2]) forming a peak in 
the surface density (i.e. a ring) in the outer disk, 
which further increases the gravitational attraction 
and results in a self-amplification of the ring for- 
mation process. Due to the growing angular mo- 
mentum of the infalling gas, the radius of the ring 
wanders outwards from i? ~ 1 kpc to i? ~ 2 kpc 
at the time of fragmentation (i ~ 100 Myr). The 
top row of Fig. m displays the surface density and 
the gravitational potential as a function of radius 
at a time shortly before fragmentation. The domi- 
nant ring corresponds to a local minimum in the 
gas-dominated gravitational potential and there- 
fore the ring is able to attract material from both 
the inner and outer regions. 

• In Model B a disk with a radius of i?disk ~ 5 kpc 
forms by simultaneous accretion of material at all 
radii R < i?disk leading to an initially fiat surface 
density distribution with an edge, with the disk 
being less distinct than in Model A. The edge acts 
again as a focal point and leads to large gravita- 
tional forces attracting material to this region and 
to a dip in the gravitational potential resulting in 
a ring at i? '^ 4.0 kpc as depicted in the middle 
row of Fig. m The sizes of the initial disk and 
the ring are determined by the angular momentum 
profile of the initial constant density sphere and the 
corresponding radii of centrifugal equilibrium. 

• In Model C the disk forms inside-out with accretion 
during the first phase occurring at a growing range 
of radii. The gravitational potential of the gas com- 
ponent is monotonically increasing at all times all 
the way from the center to the outer regions. In- 
falling material is attracted to the central regions, 
which are growing from ongoing vertical accretion 
and gas migrating radially inward. The first phase 
of gravitational fragmentation occurs in the central 
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Fig. 5. — The rotation curves for the ring systems of Models 
A,B,C at the times corresponding to Fig. [4] The vertical dashed 
lines mark the radial positions of the rings. 



region at i '^ 290 Myr, and unlike Models A and B 
not in a ring structure. In the next phase the influx 
of new material peaks at larger radii again. The 
recently accreted gas forms a ring-like structure 
around the already fragmented, irregular, massive 
gas clumps in the disk within the ring, which dom- 
inate the gravitational potential and thus influence 
the morphology of this less distinct structure (cf. 
bottom row of Fig. |4|). 

Fig. [5] shows the measured gas rotation curves for the 
three simulations at the corresponding time depicted in 
Fig. 2] with the vertical dashed lines indicating the ra- 
dial positions of the rings. Model A shows a monotoni- 
cally increasing rotation curve with a shallower increase 
in the low surface density inner region and a steeper part 
in the ring region, as one would predict from the corre- 
sponding mass proflle. Outside R ^ 2.5 kpc the rotation 
curve drops indicating that no equilibrium disk has yet 
formed in this region. The rotation curve for Model B 
exhibits less monotonic behavior indicating a clumpier 
surface density structure. The curve is steeper in the 
center resulting from a small counter-rotating gas clump, 
which does not affect the dynamics of the rest of the disk. 
Moving outwards Model B shows similarly to Model A 
a transition from a flat inner region at i? ~ 3 kpc to a 
steep ring-dominated region. At the time correspond- 
ing to Fig. |4] the disk in Model C has already frag- 
mented. As a result we see strong local variations in 
the rotation velocities produced by massive clumps dom- 
inat ing the gravitational field in their surroundings (see 
also llmmeli et al]|2004bt ). Because of the higher central 
surface densities model C displays the largest rotation 
velocities and the steepest increase in the rotation curve 
of all the three models. 

The inclusion of star formation and supernova feedback 
in Model SF adds an additional source of pressure and a 
collisionless component to the simulation. The effects of 
feedback prevent the formation of an initial distinct ring 
as seen in Model B, which has identical initial conditions. 
However, as can be seen in Fig. [3]also Model SF develops 
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ring structures, which are short-hved and less distinct 
than in Models A and B and more similar to the ring in 
Model C. 

3.3. Clumps and bulges 

Fragmentation due to gravitational instability leads to 
a clump-dominated evolution phase of our systems. The 
Toomre parameter 



Q = 



ttGS' 



(3) 



predicts instability for Q < 1, where k is the epicycle 
frequency, an is the radial velocit y dispersion and E is 
the surface density, (lToomrelll964 for a re vised criterion 
for turbulent disks see lRomeo et al.ll2010l ). 

As the temperature of the gas in the disk is typi- 
cally T ~ 10'* K, the corresponding sound speed is 
c, ~ 10 kms""'^. As we show in section [5^ the veloc- 



algorithm to detect clumps, setting a lower mass limit 
of Mniin — IO'^'^Mq, which corresponds to N^nin = 158 
SPH particles. The smoothing kernel is calculated us- 
ing a number of iVspH = 40 particles and the clumps 
are thus well resolved. We also perfor med tests using 
the clu mp detection method described in lBournaud et al] 
(|2007[ ). which defines a clump as a region with a mass 
above the minimum mass M^i^ and an over-density of 
AS > 2S(i?), where S(i?) is the average surface density 
at the disk radius R and AS = S — S(i?). Using this 
clump definition does not alter our conclusions. 

The clumps forming in the fragmentation phase have 
typical initial masses of M ^ 10^ Mp,. The Jeans mass. 



Mi 



2.92crV(GiV°'^) 



ity dispersion of the gas is significantly higher than Cg in 
the region undergoing fragmentation, which justifies the 
use of api for the calculation of Q. The velocity disper- 
sion at this phase is nearly constant with surface density 
(see section [nH]) . Thus the surface density distribution 
S, which evolves differently as a function of radius in 
the three models, determines the different fragmentation 
patterns. Consequently, fragmentation first occurs in the 
regions of highest surface density, which are the ring re- 
gions in Models A and B, and the center in Model C. 

In the left panel of Fig. IHlwe plot a map of the Toomre 
Q parameter in the disk of Model A in the ring phase at 
i '^ 95 Myr. The ring is beautifully visible as the region 
with the lowest values of Q, with Q < 1 indicating that 
it is unstable. The ring has already started to fragment, 
as the shape of the Q < 1 region is no longer circular. At 
larger radii some Q < 1 regions are also visible, however 
they are outside the equilibrium disk region and the de- 
termination of Q is thus not significant, as the criterion 
holds only for a thin, differentially rotating disk. 

After the fragmentation phase the evolution of the sys- 
tem is dominated by clumps. We use a friends-of-friends 



(jBinnev fc Tremaind [1987), 
gives an upper mass limit of masses stabilized against 
gravitational collapse by pressure support. Using the re- 
lations S w 2rhoh and h w cr^/(7rGS), where h is the 
scale-height of the disk, we find 



Mj = 1.0 X 1O*M0 



f ^- V 

V30kms-i/ 



S 



lO^Mo/kpc^ 



(4) 

Fig. |4] shows that the typical surface density in rings 
is lO^M0/kpc^. In section |331 we show that the cor- 
responding vertical velocity dispersion az in the ring 
phase is 30-40 kms""'^. Consequently, clumps with masses 
of IO^Mq are consistent with local Jeans theory. The 
clumps quickly grow in mass due to merging and accre- 
tion of diffuse material. The fraction of gas within the 
disk region, which is bound in clumps also increases with 
time reaching very high values of /dump ^ 80% for the 
gas-only simulations. 

The top left panel of Fig. [7] depicts the surface density 
map for Model A at i ~ 280 Myr, which shows four mas- 
sive clumps (M '^ 10^-^AIq) and several smaller clumps 
with masses of M ~ 10^ Mq. Two of the massive clumps 
(in the lower half of the plot) are closely interacting. In 
the right panel of Fig. [6] we plot a map of the Toomre 
Q parameter at this time. The clumps are visible as re- 
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of i = 60°; Bottom Right: Line-of sight velocity dispersion ai^^ map for an inclination angle of i = 60°. 



gions of Q < 1 . The initially rather spherical clumps de- 
velop during the simulation into rotationally-flattened, 
centrifugally supported mini-disks that interact in the 
disk plane. The global Toomre Q parameter is thus no 
longer an appropriate measure for the stability of the 
clump regions. New clumps primarily form at the outer 
radii of the disk due to newly accreted material. These 
outer regions are also Toomre unstable with Q < 1, as 
can be seen in the right panel of Fig. [51 

The mini-disks are typically aligned in the mid-plane 
with the exception of Model B, where the fragmentation 
occurs during a phase of ongoing strong accretion. The 
vertical center-of-mass velocity of these clumps is thus 
higher, which results in vertical offsets in the position 
of the clumps. However, also in Model B, the vertical 
velocity dispersion within the clumps and their vertical 
extent are as small as in Models A and C. 

The gas-only models do not show a clear bulge forma- 
tion process. This is partly caused by the fact that we did 
not run these simulations for longer than t ^ 450 Myr, 
when massive clumps have formed and we thus expect 
star formation to play a major role. In all these simula- 
tions however the most massive structure settles in the 
most central position and can thus be interpreted as a 
proto-bulge structure. In Models A and C the final stages 



show typically ^ 5 clumps with masses Afciump ^ IO^Mq 
and ~ 10 smaller clumps. Model B shows a high clump- 
merger rate in the final '~ 50 Myr resulting in only 3 
massive and three smaller clumps. Our gas-only models 
are thus different from previous models of gravitationally 
unstable, clump-dominated disks, which included pre- 
scriptions for star fo rmation (e.g. In imeli et a l. 2004^ 
IBournaud et al.|[2Q07D . These models typically show that 
massive clumps migrate inwards and coalesce to form a 
bulge within a few 100 Myr. The main reason for the 
difference is that almost all of the gas is bound in mas- 
sive clumps and the transfer of angular momentum can 
only occur between these clumps or between a clump 
and the dark matter component. The dynamical friction 
timescale is shorter for more massive clumps and thus 
they migrate to central positions more quickly. How- 
ever, our simulations distinctly lack a diffuse gas compo- 
nent which could absorb angular momentum. We have 
run isolated un s table disk galaxies similar to those of 
IBournaud et al.l ()2007f ) and were able to reproduce the 
bulge formation process by including star formation. 

As discussed in section 13.21 the effects of feedback sig- 
nificantly alter the formation of over-dense structures 
within the disk region of Model SF. The clumps in this 
model are consequently less dense, less distinct and less 
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Fig. 8. — The global radial velocity dispersion (Jji^gioh (black lines) and global vertical velocity dispersion cr^^gioh (purple lines) vs time t 
for the disk-gas in Models A,B and C. For Model A disk formation starts instantaneously, whereas for Models B and C it is delayed until 
t ~ 200 Myr and t ~ 150 Myr respectively. 

stant value of r^ 100 Mqyi'^ between t ~ 250 Myr and 
t ~ 400 Myr, which is the phase dominated by the ir- 
regular disk structures. These SFRs are thus similar to 
observations of high- 2: disks. The SFR then steadily de- 
creases due to gas consumption, being ~ 60 M0yr~^ 
at t - 500 Myr , - 30 Afoyr-^ at t - 1 Gyr and 
~ 20 M0yr-i at t - 1.5 Gyr. 



flattened than in Models A,B,C. However, as can be 
seen in Fig. |2] also Model SF develops irregular, clumpy 
morphology and ring structures. The density contrast 



AStot/Stot, where Stot = S 



gas 



of these struc- 



tures is significantly lower than in the gas-only simu- 
lations. Howev e r, if we lower the clump criterion of 
iBournaud eFall (I2007D to AS > S(i?) we stiU detect 
clumps of masses Mdump ^ lO^~^Af0. The second row 
in Fig. [3] clearly depicts two of these structures in the 
face-on stellar surface density maps. 

The morphology of the developing disk changes rapidly 
on timescales of a few t ^ 10 Myr. The d isplayed mor- 
phologi es show similarities to the models of llmmeli et al.l 
(|2004af ). The stellar component and the larger amount 
of diffuse gas in this model, allow a more efficient out- 
ward transport of angular momentum and lead to a more 
efficient inward migration of clumpy material than in 
the gas-only simul ations, in broad agr e ement with pre- 
vious work (e.g. Ilmmeli et al.l l2004at IBournaud et al.l 
l2007f ). The system of clumps initially displays various ir- 
regular morphologies before forming a bar-like structure 
{t r-^ 500 Myr). The morphology of the bar varies, but 
it persists until the end of the simulation at t = 1.5 Gyr 
as is depicted in the lower rows of Fig. [3l At the fi- 
nal stage the model does however not depict a central 
excess of surface density E, as had be en found in sim- 
ulatio ns of isolated, unstable disks by IBournaud et al.l 
(|2007( ). The radial surface density profile is well-fit by 
an exponential with a scale-length Rd ~ 2.8 kpc out 
to -R = 12 kpc. The final stellar mass of this disk is 
-M^steiiar ~ 6 X lO^'^ Mq with a gas fraction of ~ 25 per 
cent. Initial test runs with static dark potentials and 
of isolated disks produced stronger central components 
(bulges). This indicates that the interaction of the bary- 
onic component with the dark halo plays a strong role 
in reducing the centra l dark matter density (see also 
iJohansson et al.l l2009b( ) . Ongoing accretion of material 
with increasing angular momentum onto the disk can also 
prevent a central surface density excess. 

Star formation occurs mainly in the over-dense struc- 
tures. The morphologies are very similar for both the 
gas and stellar components. The SFR averaged over 25 
Myr periods, increases initially from ^ 30 MQyr"^ when 
the disk starts forming at t ^ 220 Myr to a nearly con- 



3.4. Velocity dispersion 

We define a disk region by determining a disk radius 
i?disk and a vertical extent hz selecting the rotationally 
supported part of the centrally forming baryonic object. 
We then calculate the vertical and radial velocity disper- 
sions CT_R„giob and fTz,giob asQ 

ai =< Vi>-< Vi >^ (5) 

taking into account all gas particles in this region. These 
quantities are thus measures for the global non-circular 
velocities and disk thickness. Radial inflow velocities, 
which vary with position could contribute to crK,giob- 
However, we do not observe regular inflows in our gas- 
only Models, which are dominated by clump-clump in- 
teractions after fragmentation. We plot the dispersions 
as a function of time in Fig. |S] for the Models A,B and 
C. As the temperature of the gas in the disk is typically 
T ^ 10*K, the sound speed is Cg ~ 10 kms~^. As Cg < ct 
is valid in the disk regions, we focus on the discussion of 
velocity dispersions. 

For Model A (left panel) we find that during the 
early formation phase the vertical dispersion is high 
(o'2,giob ^ 50 kms~^) and the radial dispersion stays 
constant (cr /j,giob ~ 40 kms^^) until the fragmentation 
of the ring. The vertical dispersion, a-z^gioh starts to 
decline after i ■^ 50 Myr to a nearly constant level of 
o'2,giob ^ 15 kms~^ after t ^ 300 Myr. The radial disper- 
sion, (Tfl,giob, on the other hand starts to increase with 
the onset of gravitational instability reaching values of 
ffl.giob '^80 — 90 kms^^ with strong fluctuations caused 
by major clump-clump interactions, which hardly effect 



Cz^glob- 



We ran a test-simulation of Model A without 



gas self-gravity and found that for this run the z and 

^ We use a subscript glob to indicate that the dispersion is global, 
i.e. calculated for the whole disk region. 
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R dispersions show similar declining behavior, thus con- 
firming that disk self-gravity is driving the high radial 
dispersions. 

For Model B (middle panel of Fig. [8|) the initial 
amplitude of the velocity dispersion is larger (cTgiob ~ 
80 kms~^), which is a result of the almost simultane- 
ous infall of disk material at a large range or radii. Af- 
ter this first formation phase, (Jz,gioh decreases steeply 
to reach typical values of az^gioh '^ 35 kms~^ at i ~ 
400 Myr, whereas (Ju^gioh increases to values of cr/j^giob "^ 
100— 120 kms""'^, thus resembling the evolution in Model 
A. After showing a clear decreasing tendency, the evo- 
lution of cr^^giob is affected by a temporary increase at 
t ^ 390 Myr. This increase coincides with two close in- 
teractions of two clumps each, which have an offset in 
vertical position. This in turn results in vertical gravita- 
tional forces increasing the velocity dispersion temporar- 
ily. The merger activity remains high at the end of this 
simulation preventing a further decline in cr^^giob- 

Model C features a longer phase of significant (vertical) 
accretion onto the inner disk regions and an inside-out 
growth of the disk resulting in a delay of the fragmenta- 
tion phase compared to Model B. The right panel of Fig. 
[8] depicting model C thus shows a combination of the 
features of the other two Models: cr_R,,giob exhibits first a 
relatively constant phase with ai^.giob ^ 40 kms~^, fol- 
lowed by a steep increase after the gravitational instabil- 
ity sets in. The vertical dispersion a^^gioh is initially at a 
high constant level (ctz, giob ~ 60 kms~^) before declining 
steeply to reach a typical value of az,gioh '^ 20 kms^^ in 
a similar way as in Model A. 

For the corresponding rotation-to dispersion ratios 
V^/ CTgiob, we find for radial velocity dispersion that 2 < 
V/<^R,gioh < 6 holds at all times in Mode ls A, B and 
C, wh ich is in agreement with observations (jCresci et al.l 
|2009| ). For the vertical velocity dispersion, T^/CTz_giob is 
only within this interval during the accretion dominated 
phase. After that, V/az,gioh increases to values of ~ 20, 
more similar to low-redshift disks. 

The decline in the vertical gas velocity dispersion 
o'z,g\oh displayed by Models A, B and C results from the 
dissipation of e nergy in random moti o ns E rr- As has 
been argued bv lElmegreen fc BurkertI ()2010[ ). the time 
evolution of E„ is determined by the energy input by 
accretion (or any other driver of dispersion) and the dis- 
sipation of energy, which acts on a timescale r ^ h/a, 
where h is the vertical extent of the disk. For ft, ^ 1 kpc 
and CT ~ 50 kms^^, the timescale is t '^ 20 Myr. The 
decrease in a in our models is slower, as accretion pro- 
cesses are still active, but becom e increasingly less im- 
portan t compared to dissipation. lElmegreen fc BurkertI 
(|2010f l argue, that accretion can only drive high veloc- 
ity dispersions in an initial phase of disk formation, in 
agreement with our models. 

The analysis above suggests that there should be differ- 
ent velocity dispersion patterns for face-on and inclined 
observations of our objects in the massive clump phase. 
We show in Fig. [7] the predicted observational patterns 
of Model A at t '^ 280 Myr, a time at which the ex- 
pected difference due to the inclination should be large. 
We choose Model A as the accretion in the later phases 
occurs primarily at the outer edge of the disk enabling us 
to distinguish between accretion and disk-gravity driven 
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Fig. 9. — Vertical velocity dispersion a^ vs Surface density S. 
Upper Panel: Model A at 95 Myr (black), 280 Myr (blue) and 
400 Myr (red). Middle Panel: Model B at 235 Myr (black), 375 
Myr (red) and 400 Myr (blue). Lower Panel: Model C at 250 Myr 
(black), 355 Myr (blue) and 460 Myr (black) 

velocity dispersion. Furthermore, in gas-only models we 
do not have the additional effects of star formation and 
feedback complicating the picture. The top left panel of 
Fig. [7] shows that the system is dominated by four mas- 
sive clump complexes, of which the lower two clumps are 
interacting. The presence of these clumps also affects 
the predicted line-of sight velocity map for an i = 60° 
inclined view as shown in the bottom left panel. The top 
right panel depicts the line-of-sight velocity dispersion 
map for a face-on and the right bottom for an i = 60° 
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inclined view. Thus, the top right panel shows the ver- 
tical z-dispersion. High gas surface densities correspond 
to minima in o^ ■ The impact of accretion is reflected by 
a ring-hke structure of high a^. at large radii. However, 
on the other hand, for the inclined view the high surface 
densities correspond to high values of line-of-sight disper- 
sion, due to the high radial dispersion. In contrast to the 
global quantities displayed in figure 13 these velocity dis- 
persions are local quantities, reflecting the substructure 
in the disks. 

We note that not all clumps correspond to high ur 
at all times as there are several effects playing a role. 
As discussed in section [231 the clumps are centrifugally- 
supported mini-disks showing a rotational velocity i^rot 
increasing with radius in the center-of-mass frame. This 
variation in ?;,ot contributes to the line-of sight dispersion 
uios- However, the clumps are hardly ever axisymmetric 
(cf. Fig. [71), but shaped by clump-clump interactions 
and tidal forces, which in each clump create an intrinsic 
radial velocity dispersion. Major interactions of clumps 
such as mergers produce high intrinsic planar dispersions 
thus resulting in the highest values of the inclined-view 
dispersion. A third effect contributing to the high values 
of cT/j^giob in Fig. [His the clump-to-clump dispersion orig- 
inating from the different radial center-of-mass velocities 
of the individual clumps. 

Unlike at i ~ 280 Myr (Fig. [71), where high surface 
density S is correlated with low vertical velocity disper- 
sion CTz , this is not the case for the ring phase of Model 
A at i « 95 Myr, which is displayed in the top row of 
Fig. [31 Fig. [51 shows a relatively high global CTz^giob at 
that time. In Fig. [31 we plot disk surface density E vs 
local vertical velocity dispersion cr^. The dispersion was 
determined for a pixelated map of the disk for regions 
within the disk radius i?disk for Models A,B and C at 
three different stages of evolution. The black curve in 
the upper panel depicts the ring phase of Model A, when 
dispersion is mainly driven by accretion and g^ is rel- 
atively independent of E as Cz > 40 kms~^ applies to 
a wide range of surface densities. In the regions with 
the highest surface densities, which correspond to the 
clumps that are about to form, the dispersion is already 
slightly smaller. For the subsequent clump-dominated 
phase, represented by the blue and the red curves, there 
is a clear trend of a^ decreasing with increasing surface 
density with a^, < 20 kms^^ for the high density regions, 
a trend that is also depicted in Fig. [71 

For Model B, shown in the middle panel, the initial 
disk/ring phase represented by the black curve corre- 
sponds to high vertical dispersions for all surface den- 
sities. The red and the blue curves again illustrate gen- 
eral trends of decreasing g^ with time and of lower dis- 
persion for higher densities, similar to Model A. How- 
ever these curves also show features contradicting these 
trends, most prominently the peak of the blue curve at 
t ~ 400 Myr for high surface densities. This distinct 
feature originates from the clump-clump mergers with 
vertical offset as discussed in Section 13.31 The corre- 
sponding curves for Model C are depicted in the lower 
panel. They reveal a similar evolution to those of Model 
A. 

The low 2;-dispersions in the massive clumps result 
from the efficient dissipation of energy in random mo- 
tions E„ in the dense gas. It is easier to drive and sustain 
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Fig. 10. — Intrinsic velocity dispersion (Tintrinsic vs sin(i), where 
i is the incUnation angle, under which the disk was observed. Blue 
dots are for SINS galaxies (Forster Schreiber et al. 2009i) and red 
squares are for other disks (for references see text). 

high vertical velocity dispersion by accretion in low sur- 
face density regions than in high E regions. The reason 



for this is, that E^ 



depends strongly on surface 



density, whereas the energy in the accretion flow per area 
is independent of surface density. Moreover, the accre- 
tion rates on average decline in all models resulting in a 
decline in g^ with time also for low-to medium surface 
densities. 

An anisotropy in velocity dispersions with gr, > Gz 
for gravitat i onally unstable disks was already discussed 
bv iToomrd (|1964[ ). iBottemal ([2003') found gr > g^ for 
the gas component in SPH simulations of (marginally) 
Q-stable, isolated late-type disk galaxies, but attributed 
this finding t o supern ova feedback effects. More recently, 
lAgertz et al.l ()2009aD found a similar anisotropy for HI 
gas in their study of large scale galactic turbulence in 
simulations of isolated disk galaxies. They note that 
their velocity dispersions are 'gravity driven', which in- 
cludes gravitational instability with Q < 1, as in our sim- 
ulations, but also non-axisymmetric perturbations with 
Q>1. 

In Fig. [TUl we plot observed intrinsic velocity dis- 



persions fJintr 



m massive z 



2 disk galaxies from 
the SINS survey (Forster Schreiber et al . 2009) (blue 
dots) and other hi gh-z d is ks (red squares: [Epinat et al 
20091: iStark et al.l l200l Ivan Starkenburg et all 1200 ' 



Wright et al.ll2007[ ) as a function of the inclination angle 
i, under which the disk has been observed. The obser- 
vations reveal a systematic variation of velocity disper- 
sion with sin(i), where disks that have been observed un- 
der the lowest inclination angle (sin(i) — 0.4), show the 
lowest dispersions crintrinsic- These disks typically have 
a dispersion of txintrinsic ~ 40 kms~^, whereas strongly 
inclined systems show values as high as iJintrinsic ^ 
90 knis"^. This is similar to our results displayed in 
Figs. [TlandlHl where the vertical dispersion Gz,gioh, which 
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Fig. 12. — Stellar exponential scale heights h^ vs disk radius R 
for Model SF at times t = 250, 500,800,1100,1500 Myr. 

is equal to the line-of-sight dispersion for face-on disks 
(sin(z) — 0.0) is systematically lower than the radial 
dispersion CT/j_giob, which contributes to CTintrinsic more 
strongly for more inclined disks. Unfortunately, no data 
for face-on disks is yet available, and the line-of-sight ve- 
locities of disks with sin(z) = 0.4 still have significant 
contributions from radial velocities. The lowest values 
for CTintrinsic tend to bc slightly higher than the values for 
Cz,giob in our simulations, with Model B being the best 
match to these observations. 

3.5. Stellar Disk Scale Heights 

Finally, we study whether small vertical gas velocity 
dispersions also lead to cold stellar disks with small scale 
heights. The global gas velocity dispersion CTgas,giob in 
Model SF, which is displayed in Fig. [11] initially ex- 
hibits similar behavior to Model B. However, because of 
the lack of distinct clumps and clump interactions the in- 
crease in the global radial velocity dispersion CTjj,gas,giob 
is delayed and limited to values of ^ 40 — 70 kms~^ . The 
vertical velocity dispersion CTz^gas.giob declines similarly 
to the gas-only models. 



We calculate global stellar velocity dispersions 
CTstars,giob at a Certain time T in Model SF including all 
stellar particles that have formed at t < T and show them 
in Fig. [TTJ Defined as such, the included stellar popula- 
tion changes from one time T to the next. The stellar ve- 
locity dispersion CTgtars.giob for Model SF are significantly 
larger than the gas velocity dispersions. The vertical 
stellar dispersion CTz^stars.giob is almost constant, increas- 
ing from ^ 60 kms""'^ at t ~ 400 Myr to ~ 75 kms~^ at 
t ~ 1.5 Gyr. CTi?_stars,giob howcver initially (t < 500 Myr) 
shows a similar trend to CTfl_gas,giob with the dispersion 
rising to values of CT/j^stars.giob ^ 90—100 kms~^ in the bar 
phase after a minimum of CT^^stars.giob ^ 60 kms~^. The 
evolution of stellar velocity dispersions CTstars.giob with 
ti me t as displayed in F ig. [TTjis thus similar to the results 
of llmmeli et al.l (|2004at ) for central stellar velocity disper- 
sions in simulations of gravitationally unstable disks. 

In Fig. [12] we display stellar exponential scale-heights 
hz as a function of the disk-radius R extracted from 
an edge-on view of the galaxy. In the initial formation 
phase the stellar disk is thick with stellar scale heights 
of hz ^ 600 — 800 pc as can also be seen in the top 
row of Fig. [3] In this phase the disk is dominated by 
strong accretion- processes as discussed in section 13.41 A 
population of stars that have formed in clouds before en- 
tering the disk settle in a thick disk. However, new stars 
typically form near the mid-plane of the disk and con- 
sequently the scale-heights decrease in the phases with 
the highest SFRs. Typical values are hz '^ 200 - 400 pc 
for t ~ 500 - 800 Myr, comparable to the local Milky 
Way thin disk. This evolution is in agreement with a 
constant CTz^stars.giob, as the surface density of the disk 
is continuously increasing. After t ^ 1 Gyr the scale- 
heights in the central region begin to increase again up 
to values hz ~ 500 pc. This heating is caused by the bar 
structure. The scale-height hz depends only weakly on 
disk-radius with the strong early variation at 250 Myr 
caused by the small radial extent of the disk, which is in 
the first stages of its formation. Finally, we note, that 
the specific prescription for star formation and the choice 
of the star formation density threshold pth, influences the 
formation of the disk, especially in the early phases, and 
thus the disk scale heights. 

4. DISCUSSION 

We have shown that different infall scenarios with disks 
forming inside-out or simultaneously at a large range of 
radii can result in the formation of rings. We argue that 
distinct edges of disks lead to enhanced accretion of ma- 
terial at these locations which results in rings. Such fo- 
cal points could also originate from enhanced infall of 
material with a certain angular momentum or scenarios 
which have less continuous infall, for examp l e an evo- 
lutio n driven by minor mergers ()Naab et al.l 120091 see 
also iBurkert fc HartmannI l2004f ) . The position of the 
ring is initially determined by the angular momentum 
distribution of the material, which forms the first disk 
in the center of the halo. As displayed by one of our 
models, the ring radius can increases later-on if material 
with higher angular momentum settles into the equato- 
rial plane before the ring breaks up. We conclude that 
these ring formation processes should play a role in mas- 
sive high-redshift ring galaxies. We also find that rings 
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and clumps characteristically influence the shape of the 
rotation curves of galaxies. 

Our simulations indicate that observations of rings, 
clumps and bulges in disk galaxies at redshift z ^ 2 are 
likely to display systems during different stages of evo- 
lution. Our gas-only models are not able to reproduce 
the formation of bulges very well, as they lack diffuse 
gas and stellar components. However we expect mas- 
sive star formation to set in once massive clumps have 
formed, which would lead to a faster migration of clumps 
to the center as shown in Model SF. When we include 
a simple, but widely-used model for star formation and 
subsequent feedback, our simulations still display the for- 
mation of clumpy and ring-like structure, however their 
density-contrasts compared to the rest of the disk are 
lower than in the gas-only simulations. Stars mainly 
form in these structures, which are not long-lived. Our 
treatment is too simplistic to allow for a statement on 
whether feedback can destroy clumps in z ^ 2 disks (cf. 
IKrumholz fc Dekell[2010l: IMurrav et alllMol) . 

Interestingly, our Model SF does not show a strong cen- 
tral component in its final stage, its radial surface density 
profile is relatively well fit by an exponential with a scale 
length comparable to that of the Milky Way. Transfer of 
angular momentum from the baryonic component to the 
dark matter plays an ess ential role, whi c h is in agreement 
with the conclusions of lBurkert et al.l (I2010D. who con- 
clude d that adiabatic contraction (e.g. iBlumenthal et all 
Il986f l of the halo due to baryonic infall is problematic, 
as it predicts halo spin parameters A that are signifi- 
cantly higher than ex pected. It is also in agreement with 
iPedrosa et al.l (|201Clf ). whose simulations of galaxy for- 
mation in dark halos are also not able to reproduce the 
predictions of the adiabatic contraction theory. 

We have run Model SF without feedback and find that 
due to the lack of additional pressure, clumps with sim- 
ilar masses to those of the gas-only runs (Mdump ^ 
IQ^-^-^Mq), but a less flattened morphology, form in 
the early disk formation phase {t ~ 220 — 300 Myr). The 
star formation occurs mainly within them at extremely 
high rates (5'i^i?initiai ^ 300 Moyr"^) consuming almost 
all first-infall gas within ~ 150 Myr. After a clump- 
growth and migration phase the clump phase lasts until 
t ~ 450 Myr. The final system again does not display 
a strong central component, the bulge is however more 
massive than in Model SF. We note that the observa- 
tional signatures of typical clumpy z ^ 2 disks (SFR, 
clump fraction) are intermediate to our simulations with 
and without feedback. 

We find that accretion of gas can drive high velocity 
dispersion of the order of 40 — 70 kms~^ in radial and 
vertical direction, as long as it penetrates disk regions of 
low to medium surface densities. Thus accretion plays 
an important role in the transient early phases of disk 
assembly, but affects the massive clumps region in later 
evolutionary stages less strongly. We note, that the in- 
fall of gas clouds produced by thermal instability onto 
the disks in our simulations enhances high velocity dis- 
persions. The thermal instability arises due to numerical 
noise in the initial conditions, h owever, as has been ar- 
gued by iKaufmann et al.1 ()2006D , realistic counterparts 
of these clouds are likely to exist in gas accretion onto 
galaxies. The decrease of cr^^giob with time is consistent 



with est imated dissipation timescale s and with predic- 
tions of lElmegreen fc Burkertl ()2010f ) Gravitational in- 
stability and the gravitational interactions between mas- 
sive clumps however are able to produce high dispersions 
for in-plane motions. The cr^, giob vs time plot of Fig. [8] 
reveals phases of driving velocity dispersion by accretion 
an d disk- instability, in agreern ent with the predictions 
of lElmegreen fc Burkertl (|2010D . Peaks of cr n^giob pro- 
duced by close interactions of the most massive clumps. 
Clump-clump mergers increase the radial dispersion for 
all our models. However, only in Model B, which fea- 
tures the largest vertical offsets of clumps, they also effect 
the vertical dispersion due to vertical gravitational forces 
between vertically offset interacting clumps. Except for 
phases dominated by mergers of this kind, all models 
show a characteristic decrease of cr^^giob with time t, as 
the dissipation of energy in random motions starts to 
become increasingly more effective and to dominate over 
the driving of high velocity dispersions by accretion. The 
rotation-to-radial gas velocity dispersion ratios V/ajigioiy 
agree very well with observations of high-z disks for our 
Models A, B and C at all times. For the vertical disper- 
sion, V/az,gioh only agrees with observations during the 
accretion-dominated phase. 

In the clump dominated phase, all models also show 
a characteristic decrease of CTz with increasing surface 
density S , as long as they are not affected by mergers 
of clumps with vertical offsets. Because of this, face- 
on views of clumpy disks display minima of the line-of- 
sight dispersion at the positions of the clumps, whereas 
for an inclined view, the clumps correspond to higher 
than average dispersion resulting from the high in-plane 
dispersion. This is not true for the ring phase, in which 
ffz is relatively independent of E. However, already in the 
phase of initial ring fragmentation, the clumps, which are 
about to form, show smaller values of CTz than the rest of 
the ring. 

We present observational evidence for a systematic 
variation of vertical velocity dispersion with viewing- 
angle of the disk for z ~ 2 galaxies. The observed 
dispersion for disks with the lowest inclination angle 
sin{i) ~ 0.4 is still higher than the vertical dispersion 
in our Models A and C, only partly explained by the 
contribution from in-plane motions. Model B is a bet- 
ter match to the observations for low-inclination disks. 
A possible explanation might be that feedback processes 
and/or ongoing strong vertical accretion into the massive 
clump region create high radial and vertical velocity dis- 
persions of the order of the observed face-on value and 
vertical offsets of clumps and thus play a significant role 
vertically, but not for the in-plane motions, whose dis- 
persion is dominated by disk self-gravity. The symmetry 
of our initial conditions also defines a perfect disk plane. 
Mini-disks inclined to the disk plane would also give rise 
to vertical dispersions. 

The scale-heights of our forming stellar disk in Model 
SF are significantly t hinner than those simulated by 
iBournaud et al.l (|2009( ) and do actually not correspond 
to thick disks. Only the early phase, which is character- 
ized by high accretion rates into all disk regions displays 
exponential scale-heights h^ ^ 600 — 800 pc. The fi- 
nal scale-heights hz ~ 500 pc are larger than that of the 
Milky Way thin disk but lower than that of the thick 
disk. One reason for this is that our Model SF includes 
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feedback and does not display distinct, massive and long- 
lived clumps. In our test-run without feedback, we find 
that scale-heights are higher by a factor of ~ 1.5 due 
to the more compact clumps that scatter stars more ef- 
ficiently, and show stronger radial variations, but still 
they are thinner than the Milky Way thick disk. We 
have also run test-simulations of isolated, massiv e disks 
with similar initial conditions to iBournaud et al.l ()2007f ) 
and have been able to produce scale-heights of the order 
of 1 kpc. The final scale-height depends critically on the 
Toomre Q factor of the initial disk. The lower Q, the 
more massive and distinct the clumps that form and the 
more vertically extended the disks. 

The scale- heights are however always smalle r by a fac- 
tor of ~ 2 compared to IBournaud et all (|2OO90 . The rea- 
son for this might be that Bournaud et al. used a sticky 
particle code for the gas component, whereas we use an 
SPH code. We also find small variations of scale-height 
with radius, possibly on account of different halo poten- 
tials. The increase in scale-height is consistent with the 
idea that the disk thickens by stars being scattered off 
massive clumps. In contrast to collisionless stellar par- 
ticles this is not possible for massive gas clouds, which 
due to their larger cross-section dissipate kinetic energy 
in close encounters. A simple estimate for the velocity a 
particle can gain perpendicular to its original direction 



is given by 



Ay = 



2GMb 



Voib^ 



(6) 



where b is the impact parameter, a is the scale-radius of 
the clump, Vq is the relative velocity of the clump and the 
par ticle and AI is the clump mass ( Binnev fc Tremaind 
1198?. p. 474). Taking this as an estimate for velocity dis- 
persion and using AI — 10^ Mq, a — 1 kpc, 6 = 3 kpc 
and Vq = 50 kms~^, we find AV ^ 50 kms~^. Thus disk 
heating is likely to occur in clumpy, massive galaxies. 
However, the efficiency of this process depends on how 
compact the clumps are, how long they survive and where 
and when the stars form. Because of this it is not clear, 
whether this is the only process at work in creating thick 
disks in high-redshift galaxies and whether they are the 
birthplace of present-day thick disks. Distinct, massive 
clumps are also present in our run of Model SF without 
feedback. However, the continuous accretion of gas con- 
tinuously allows continuous star formation and reduces 
the disk heating. 

Given the limitations inherent in our models, we are 
able to study essential physical processes that are im- 
portant for the early stages of galaxy evolution and their 
signatures. We are also able to qualitatively reproduce 
observational signatures of z ~ 2 galaxies allowing us to 
draw conclusions about the role these processes play in 
these objects. However, more detailed effort is necessary 
to understand in detail the observational properties of 
high-redshift disk galaxies. 
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